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ABSTRACT 


A finite difference scheme is developed for constructing quasi conformal mappings 
for arbitrary simply and doubly-connected regions. Computational grids are 
generated to reduce elliptic equations to canonical form. Examples of con- 
formal mappings on surfaces are also included. 
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1 . Introduction 

The motivation for constructing quasi conformal mappings lies in their 
application to the generation of curvilinear coordinate systems. Quasicon- 
fortnal mappings may be used to reduce any second order linear elliptic partial 
differential equation to canonical form (i.e., the principal part of the 

i 

differential operator reduces to the Laplacian). Consequently, if one is 
solving an elliptic boundary value problem, an appropriate quasi conformal 
mapping could be used to simultaneously fit the boundary contours with coordi- 
nate lines and simplify the original partial differential equation. The equation, 
in canonical form, could possibly be solved more efficiently or by methods 
which would not be applicable to the original equation in cartesian coordinates. 

A related application of quasi conformal mappings is in the construction of con- 
formal mappings on surfaces. Although they have been studied extensively by 
complex analysis, very little work has been done on the numerical construction 
of quasi conformal mappings. Of the methods which have been proposed, that of 
Bel inski i et al. [2] uses a fixed boundary correspondence which determines 
the mapping parameters, and the method of Mas ti ri and Thompson [4] would be dif- 
ficult to implement on arbitrary regions. A finite element version of the 
latter method developed by Weisel [9] appears promising. However, the class of 
mappings and the type of regions presented in the examples are very limited. 

In recent years a finite-difference method for constructing conformal 
mappings developed by Allen [1] has been used by Mobley and Stewart [5], Pope [6], 
and Yen and Lee [8] in the construction of orthogonal coordinate systems. 

Although all of these authors use essentially the same numerical method, there 
are differences in the way the boundary values and the conformal module of the 
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region are computed. This method is not as accurate or efficient as other 
conformal mapping methods using integral equations or series expansions, but 
it does have the advantage of simplicity since the module of the region, the 
boundary correspondence, and the interior grid points are determined in 
a single iterative procedure. Modeling this conformal mapping procedure, it 
will be shown how quasiconformal mappings can be constructed and applied to 
the reduction of elliptic equations to canonical form and the construction 
of conformal mappings on surfaces. Except for the method of Godunov and 
Prokopoy [3], this appears to be the only conformal mapping method which can 
be easily adapted to handle the problem of constructing quasiconformal map- 
pings, 

2. Boundary Value Problem for Quasiconformal Mappings 

Let D be a bounded simply-connected region in the xy-plane whose boundary 
C is a simple closed contour. Let Zy z 2 , z 3 and z^ be distinct boundary 
points ordered by the orientation on C. There exists a unique quasiconformal 
mapping of D onto the interior of a rectangle such that the points z- map to 
the vertices which are also ordered by the orientation of the rectangle. The 
ratio of the length to the width of the rectangle is a quasiconformal invariant 
called the module of D and will be denoted by m. The quasiconformal mapping of 
the region D onto a rectangular region can be obtained by constructing the 
mapping of D onto a square region S which satisfies the linear system 


? x = m ( CT1 y + H) 


~ -ni(an x + bn y ) 
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where ac - b = 1 . On setting (y,v) - (s,mn)> it is obvious that y and v, 
as functions of x and y, satisfy the Beltrami system and hence we arrive at 
the desired quasi conformal mapping (see [4] for further details). 

It is easily shQwn that x and y, as functions of s and n, satisfy 


ax - 23X + yx - J 2 (a + b ) 

££ Sn nn x y 

ay - 2ey + vy - J 2 (b + c ) (2) 

Sn T nn v x y J ' 7 


where 


2 2 
a = ay - 2bx y + cx 
n n n n 


S = ay 5 y n - Mx^ + x^) + cx^ 
Y = ay^ - 2bx 5 y { + cx* 


J = V* - Vr 


At each boundary point of D one of the functions, % or n, is constant 
while the other satisfies an oblique derivative condition. This implies that 
the condition 3=0 must be satisfied by x and y at each boundary point of 
S (except for vertices). Note that in the solution of (1) we would have 
3=0 throughout S. Thus the computed value of 3 can serve as a test of the 
accuracy of our solution. It also follows from (1) that m = A/y. Therefore, 
the equations (2) can be written as 


m 2 x 

2 

m y 


SS 

U 


+ x 
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= mJ(a x + b y ) 

= mJ(b x + c y ). 


(3) 
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The same procedure is used for doubly-connected regions with a periodicity 
condition applied on two opposite sides of S. 

The system (2) or (3) may be solved using an iterative procedure with 

(1) the right-hand sides of the equations, (ii) the boundary values of x and 
y, and (iii) either a, 3 , y, or an approximation of m, re-evaluated at each 
iteration. As in the previously discussed conformal mapping methods, both 

(2) and (3) have performed equally well in numerical examples. 

3 . Reduction of E lliptic Equations to Canonical Form 

The application of quasi conformal mappings in the solution of elliptic 
equations is well-known, An elliptic equation of the form 

2 

au + 2bu + cu - f (u , u ,u w ) , ac - b =1 
xx xy yy ' x y 


transforms to an equation of the form 


m 2 u 




under the transformation defined by (1). This transformation- will make'the 
solution of many problems much easier. For example, if a, b, and c are con- 
stants and f = 0, then it can be shown that g = 0 and we need only solve 

m^u „ + u = 0 on a square region. This is efficiently done by separation 
55 nn 

of variables or a direct numerical method. Computational grids for solving 
this problem on a simply and doubly-connected region are given in Figure 1. 

The practicality of using quasi conformal mappings in solving elliptic 
equations will be further examined in the following example. The function 
u = cos(x - y) is a solution of the partial differential equation 

u + u + u + cos(x - y) = 0. 
xx xy yy 


( 4 ) 


We will solve this equation numerically for 0 <_ x, y <_ u with Dirichet boundary 
conditions prescribed by the known solution. In terms of curvilinear coordi- 
nates, this equation can also be written as 


m 2 tl , ,, . 2Jm 

m u ~ + u + — 

nn yg- 


cos(x - y) = 0 , 


The quasi conformal grid for solving (5) is illustrated in Figure 2. For the 
purpose of assessing the influence of the error in the iterative solution of 
(2) on the error in the solution of (5), the iteration was stopped occasionally 
and the solution of (5) was computed. A comparison of a normalized value of 
3 with the error in the solution of (5) is plotted in Figure 3. Note that 
3 has been normalized so that for the construction of conformal mappings 
(a = c = 1 , b = 0 ), the values, along the abscissa would represent the degree of 
nonorthogonality. Equation (4) was also solved on a uniform rectangular mesh 
with the same number of grid points. The error in this solution serves as an 
approximation of the discretization error which would result in solving (5) with 
the exact quasiconformal mapping. It also serves as a test of our method against 
the traditional method for solving (4). Figure 3 indicates a nearly linear re- 
lation between the plotted variables. This is to be expected since the major 
part of the truncation error for larger 3 , is due to the omission of the mixed 
derivative term, which is a linear function of 3 . 

A few remarks concerning the numerical solution of ( 2 ) and (5) will be made. 
The system (2) was solved using point SOR in the same way one would construct 
a conformal mapping (see [7]). For the condition 3 = 0 on the boundary, a 
form of one-side "upwind" differencing was necessary to maintain convergence 
for the value of b = 1//T in this example. The elliptic equation (5) was solved 
using a direct elliptic solver (see [7]). After 35 iterations of (2) the 



maximum error in the solution of (5) was within 125 per cent of what we 
estimated to be the maximum discretization error. At this point the value of 
| 3 1 was still decreasing, but at a very slow rate. The exact value for m 
in this example is 1 due to symmetry. The computed estimate, which was 
the root-mean-square value of /a77, was 1.00006 after 35 iterations. 

This example does not illustrate an efficient use of quasiconformal 
mappings. The method would be very efficient when one had to solve an 
elliptic equation with many different boundary conditions or inhomogeous 
terms. In that case the quasiconformal mapping would only have to be con- 
structed once. 

4. Conformal Mapping on Surfaces 

A second area of application of quasiconformal mappings is in the con- 
struction of conformal (isothermal) coordinates on a surface. Let M be 
a smooth bounded surface in xyz - space which is defined by the parametric 
equations 

x = x(cj>,e), y = y(M), z = 

The parameter region "n the <|>e - plane may be an arbitrarily shaped simply or 
doubly-connected region but it is assumed that the boundary is composed of 
simple closed contours and the mapping from the parameter region to the sur- 
face has a nonvanishing Jacobian. 

A conformal mapping of M onto a rectangular region can be constructed by 
constructing a mapping from a square region of the 5n - plane onto M which 
satifies 
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where P ** (x,y,z) and m is the module of the surface M. If these equations 
are written in terms of the parametric variables ^ and 0 , we conclude that 
<j> and e satisfy 


me,. = bs - c<}< 

5 Ti n 


md)^ = ae - bf 
5 n n 


where 


a - |P 6 l Z /d 

b % ' V d 

c - |P + | 2 /d 



However, this is equivalent to (1) with (x,y) replaced by ( , 0 ) . In this case 
the quantity which corresponds to B/voy would be the cosine of the angle be- 
tween a 5 = constant and a n = constant coordinate line on the surface M. 

Conformal grids have been constructed for several simply connected sur- 
faces. Three surfaces are listed below. In the first two cases, the parametric 
region was the projection of the surface onto the xy - plane. 

i. Paraboloid: z = 1 - x 2 - y 2 , x 2 + y 2 £ 1 

2 3 

i i • Bicubic: z = xy , - 1 < x,y < 1 

iii. Torus: x = (2 + sind>) cos e 

y = 1 - ( 2^/tt) 2 . 
z = (2 + sin<}>) sin e, 

- tt/ 2 < Cj) <_ tt/2, 0 £ 0 <_ TT 
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The plots of these surfaces appear in Figure 4. It is difficult to visualize 
the orthogonality from the plots, but the departure from orthogonality was 
less than one degree except near vertices on the boundary where the ortho- 
gonality condition was not imposed. 

The advantages of conformal coordinates are well-known. Problems in- 
volving heat conduction, ideal fluid flow, and electric fields can be solved 
as easily on the surface as they can on a rectangular region in the cartesian 
plane. 

6 , Conclusions and Discussion 

A finite difference method, which has been widely used for the construction 
of conformal mappings, has been generalized to construct quasiconformal mappings,. 
This development will increase the class of problems which can be solved using 
the currently available fast elliptic solvers developed by Swarztrauber and 
Sweet [7], Even when iterative methods are required, the absence of a mixed 
derivative and a rectangular region both would tend to give faster convergence 
especially when optimal iteration parameters are known. 

We will conclude this report with an open problem. It is known that if 
a = c = 1 and b = 0, then the solution of (2) which satisfies 3 = 0 on the 
boundary of the square S will also satisfy the system (1), and hence determines 
a conformal mapping. This follows directly once it is noted that the quotient 

"v + lT| x 

is an analytic function. Here it has been assumed, and numerical results tend 
to verify, that the same result holds for arbitrary quasiconformal mappings. 
However, no proof has been found. 
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Figure 1. Quasi conformal grids for simply and doubly-connected regions, 
a * 1, b » 1/2, c - 5/4. 



Figure 2. Quasi conformal grid for the solution of elliptic equation, 
a = c « 2/v'3, b = 1//3. 







